Decoupled Spin Chains in Ba3CuSb209 
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Inspired by recent measurements on Ba3CuSb2 09, we explore the interplay of spin and orbital 
degrees of freedom in a triangular lattice for possible spin liquid-like behaviour. Our effective Ifamil- 
tonian describes the magnetic Cu layer by considering superexchange paths along the intermediate 
O's. A mean field study of ordered states reveals an orbitally ordered ground state that isolates the 
spin degrees of freedom into spin chains. This scenario naturally explains the linear specific heat as 
T — > and, by including the effects of additional defect Cu's sitting off the triangular planes, we 
find remarkable agreement with the measured magnetic susceptibility. 
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The quantum spin liquid is an exciting and heavily 
sought ground state in which frustrated interactions lead 
to novel quantum order. In ID, the antiferromagnetic 
(afm) spin chain has a spin liquid ground state and, in- 
deed, it is well understood and has been realized exper- 
imentally [9l [111 [13] ; however, the real excitement is in 
finding a gapless two dimensional quantum spin liquid. 
The prototypical frustrated system in 2d is the afm tri- 
angular lattice: with isotropic exchange, however, the 
ground state is not a spin liquid, but rather, it stabilizes 
into the 120° ordered state [5- There exist numerous 
candidate quantum spin liquids [3], including a variety 
of organic compounds and, more recently, the inorganic 
Ba3CuSb209 [T] and its sister compound Ba3NiSb209 
[2]. The magnetic ions in these compounds, the spin 1/2 
Cu and spin 1 Ni, respectively, form a triangular lattice 
in planes that are effectively isolated by two layers of 
the non-magnetic Sb. Recent measurements on the Cu 
system provide evidence that the low temperature state 
of Ba3CuSb209 may be a quantum spin liquid. Moti- 
vated by these results, we seek to explain the spin liquid 
behaviour by examining the superexchange in the trian- 
gular lattice including the two-fold orbital degeneracy of 
the Cu. 

The geometric frustration of the triangular lattice may 
be enhanced by the competition of the spin and orbital 
degrees of freedom [5J[7]; however, the coupling with the 
orbital degree of freedom often relieves the frustration 
via orbital ordering [5l[8]. Indeed, our analysis predicts a 
constant orbital ordering that entirely decouples the spin 
degrees of freedom on the triangular lattice into afm spin 
chains. ID spin chains tend to a linear specific heat at low 
temperatures and, furthermore, the magnetic susceptibil- 
ity does not show a conventional magnetic ordering (but, 
rather, has a distinctive maxima at fc^T w Jafm [S]). 
When we account for the impurity Cu substitutions into 
neighbouring layers of Sb, we naturally reproduce the 
temperature dependence of the magnetic susceptibility. 
We predict a defect substitution of Sb by 8.5% Cu (more 
than the 5% Cu substitutions reported in [1]), coupling 
/erromagnetically with the nearest Cu in the triangular 
lattice, while the in-chain Cu's couple antiferromagnet- 



ically with an exchange Jafm ~ 60 K, to be compared 
with the reported 32 K extracted for a triangular lattice 
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The Cu ions sit within an octahedron of O's, tilted rel- 
ative to the neighbouring octahedra, so that the single d 
hole lies in the degenerate e^ orbital doublet. Superex- 
change most likely occurs with an intermediate double 
occupancy of holes of an to a neutral state. Virtual 
hopping proceeds along pairs of O in nearest proximity 
of the neighbouring tilted octahedra. The hybridization 
energy of Sb is much higher so that they aren't expected 
to participate. 

Figure [T] depicts the triangular lattice of Cu's within a 
cage of O's connected to the inequivalent up versus down 
Sb layers. The Bg orbitals are defined relative to the 
(arbitrarily chosen) z-axis and will be denoted d^z^^r^ 
(or more simply just dz2), or by the spinor = 4-1, and 
dx2^y2, corresponding to = —1. The spin is denoted 
by s. For hopping along two intermediate O's directed 
along the a and b axes, the corresponding Cu-Cu 'bond' 
is labelled c ^ a,b where a,b,c G {x, y, z\. The hopping 
into the upward plane is along an ordered pair of positive 
axes, say (a, 6), while the alternative hopping connecting 
the same Cu-Cu pair into the downward plane is along 
{-b,-a). 

Assuming the full octahedral symmetry in the Cu-0 
arrangement, the intermediate hopping integrals share 
the following relations [T3]: 
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A trigonal distortion, stretching or compressing along the 
(111) direction, would not lift the e^ orbital degeneracy; 
however, these relations would become merely approxi- 
mate. Allowing for inequivalent energetics and hopping 
in the two Sb layers, we define two overall exchange pa- 
rameters Ji and J2 for exchange via the Sb(l) and Sb(2) 
layers, respectively. Since the hopping is essentially lon- 
gitudinal (the d and p orbitals are roughly aligned), we 
consider double occupancy of the same p orbital and the 
superexchange is antiferromagnetic. 

The superexchange path between layers of Cu differs 
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FIG. 1: A single Cu (solid) layer and the adjacent Sb(l) 
(striped) layer - the adjacent Sb(2) layer resides directly un- 
der the Cu triangular lattice. The O (hollow) reside midway 
between each Cu-Sb. The axes x, y and z form a right-angled 
coordinate system: the Sb(l)'s are drawn in the layer above 
the Cu layer so that the labelled local axes come diagonally 
out of the page to three such Sb(l)'s; the negative Cartesian 
directions are toward the O's sitting at right angles to the 
Sb(2) beneath each Cu. 



from the nearest neighbour intra-layer exchange path by 
additional hopping along a third O; the orientations of 
the hopping along two of the three O's is identical that 
of the hopping involved in the intra-layer superexchange. 
We can conclude that the inter-layer exchange is strictly 
smaller than the intra-layer exchange (by two orders in 
perturbation theory). As an aside, we point out that 
although the isostructural 6H-B Ni-based compound has 
no orbital degeneracy - it has been lifted by the strong 
Hund's coupling that leads to a spin 1 degree of freedom - 
the same conclusion holds for it: the intra-layer exchange 
should be strictly larger than the inter-layer exchange, in 
disagreement with a recent proposal [K] . 
The resulting effective Hamiltonian is 

^e// = ' ^i+y (Ji^i'^i^'^i+y] + J2y[Ti+y,'ri]) 

i 

+Si-S,+4(Ji+ J2)Z'^[r„T,+^] (2) 
-VSiJi^J2)Z^[n,T,+,]] 



where the orbital dependence of the various bonds is 
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Although we artificially break the full octahedral symme- 
try in our choice of e^ basis, this effective Hamiltonian 



retains the full symmetry that can be verified by a ro- 
tation of the orbital basis, such as {x,y,z) — )■ {z',x',y') 
which corresponds to a 60° rotation of the orbital spinor 
r. 

We analyze the effective Hamiltonian by assuming a 
mean field decoupling of the spin and orbital degrees of 
freedom. The spin and orbital directions are specified by 
a pair of polar winding angles aj and l3j, respectively: 



Si = s{cos{9s + a ■ Ui), 
Ti = (cos(6'^ +/3 • ai),e 
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The mean field solutions do not depend on the azimuthal 
angles, </), a', or /?': we consider rotations about the y 
axis without loss of generality. Note that the mean field 
energy never depends on the angle Og but depends on O-r 
in the special cases of 'ferromagnetic (fm)', /3 = 0, or 
'afm', /3 = TT, (or alternating) orbital ordering. 

The mean field ground state corresponds to a 'fm' or- 
bital ordering in the dx2-y2 orbital that decouples the 
spin degree of freedom into afm spin chains along the 
z axis. This can be seen immediately from the Cu-0 
arrangement depicted in Figure [TJ exchange coupling 
along the Cu-Cu x and y bonds inevitably involves hop- 
ping along the z axis along which the planar dx^-y^ or- 
bital has zero overlap. There are also degenerate solu- 
tions corresponding to chains along the x and y axes 
that are most easily constructed by a coordinate rota- 
tion in orbital space. The afm chain exchange is com- 
pletely independent of the up/down asymmetry: indeed, 
Jafm — 9(Ji -I- J2). No other mean field orbital ordering 
described by (6][7) locally minimizes the energy. 

A trigonal distortion modifies our exact relations 
amongst the hopping integrals ([!]) and hence also mod- 
ifies the precise numerical prefactors in the effective 
Hamiltonian ([2]). However, we expect our conclusions 
will remain unchanged. This is because we can rotate 
to a new basis for the e^ orbitals where the new 'pla- 
nar' orbital is chosen such that the overlap with the p 
orbital along the z-axis vanishes in the distorted octa- 
hedron. Then we are assured that an orbital ordering 
into this orbital will still break the triangular lattice into 
independent afm spin chains. 

Another possible distortion of the lattice common in 
afm spin chains is a Peierls distortion which results in 
alternating short and long Cu-Cu 'bonds' and favours 
singlet pairing of spins along the short 'bond'. This sort 
of distortion would mix the e^ orbital states and our sim- 
ple symmetry arguments break down. Without a quanti- 
fied analysis of the multiple exchange paths incorporating 
such deformations, we cannot comment on possibility of 
such a transition. 

In the samples of Ba3CuSb209, Zhou et. al. report 
that the Cu impurities are substituted specifically onto 
the Sb(2) sites [T. Such an impurity Cu is connected via 
three alternative O's (in Figure [I] these O's correspond 
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to each of the negative axes) to the nearest in-plane Cu: 
the angle formed between either Cu and the intermedi- 
ate O is a right angle (or nearly so) so that the exchange 
between the two is ferromagnetic [TT]. At low energies, 
these fm coupled Cu's will form an effective spin 1 de- 
fect. The resulting effect on the remaining afm chains 
now depends on the effective exchange of this spin 1 with 
the next-nearest Cu's on either side in the chain. Note 
that since we must consider additional exchange paths 
directly between the defect Cu and these next-nearest 
neighbour Cu's, we cannot conclude a priori if the result- 
ing exchange will be antiferro- or ferromagnetic. With 
an afm exchange, the two next-nearest neighbour Cu's 
generically screen the spin 1 defect leaving no orphan 
spin and a cut spin chain less three spin sites altogether 
[H] ; possible fine-tuning of the afm exchange leads only 
to partial screening of the spin 1 to an orphan spin 1/2 
whilst healing the spin chain so that it remains uncut 
[13j . A fm exchange with the two chain ends renormal- 
izes to an orphan spin 1, cutting the chain in two with 
one less spin. 

In order to discern our most likely defect scenario, con- 
sider the source of 1/T divergences in the magnetic sus- 
ceptibility. An orphan spin S contribution follows the 
simple Curie law: xs = [giiafSiS + l)/3kBT. The 
same density of orphan spins for spin 1 versus spin 1 /2 
therefore contributes a factor 8/3 larger l/T divergence 
in the magnetic susceptibility. Finite odd length afm 
chains also diverge as T — 0; however, the region of 
1/T behaviour is limited to ksT < 0.15Ja/m [9|. This 
difference in temperature scale allows us to separate the 
orphan spins and the broken chain contributions to the 
divergence. As a first approximation, we can extract the 
defect concentration ignoring the chain contribution: for 
orphan spin 1, we require a concentration of 8.5% defects, 
compared with ^ 20% spin 1/2 orphan spins. Although 
even this smaller estimate exceeds the reported 5% Cu 
substitutions, we will hereafter assume the spin 1 orphan 
scenario, which, recall, also entails breaking the chains 
at the site of each defect. 

The broken afm chain contribution is calculated as- 
suming a random distribution of defects. The corre- 
sponding probability that a Cu on the triangular lattice 
should find itself in a chain of length n is 

P^{n) ^nx^{l-xy' (8) 

where x is the density of defects [SJ. We fit the experi- 
mental magnetic susceptibility by combining the contri- 
butions of fm coupled pairs of spin 1 /2 defects and broken 
afm chains varying the density of defects and both the 
chain afm exchange and fm defect exchange. The finite 
chain susceptibility and specific heat was calculated by 
exact diagonalization and interpolated in 1/n for large n 
using the infinite chain Bonner-Fisher curve ^9: . 

The fit of our model with the experimental data is best 
shown by subtracting the modelled orphan spin contri- 




ksT/J 



FIG. 2: A comparison of the experimental nx (after or- 
phan 5* = 1 spin substraction) with the modelled broken afm 
spin chain response Xchains- Note that the T — > divergence 
is now due entirely to the odd length finite chains. The fit 
parameters involved a Cu defect substitution on Sb(2) sites 
of 8.5% with fm exchange with the nearest Cu in-plane of 
Jfm = 34 K, and a chain afm exchange of Jafm = 63 K. In- 
set: a graphical depiction of the triangular lattice decoupled 
into broken spin chains with spin 1 orphans at each substi- 
tuted Cu defect. 



bution and comparing Xexp — Xs=i with the broken chain 
susceptibility (see Figure [2]) fTBl . The fit parameters are 
a fm exchange of Jfm = 34 K coupling the 8.5% sub- 
stituted Cu's with the nearest in-plane Cu, and an afm 
exchange of Jafm = 63 K within the broken chains. Af- 
ter subtracting the orphan spin contribution, the T — > 
divergence is now due solely to the odd length broken 
chains. The height and position of the local maxima are 
essentially determined by the afm exchange within the 
chains. The fm exchange affects the magnetic suscep- 
tibility curve far less drastically: in practice, all three 
parameters, x, Jafm, and Jfm are fit simultaneously in a 
least squares minimization scheme. 

With these same fit parameters, we can compare our 
model's predictions against the experimental magnetic 
heat capacity. We reproduce the same linear behaviour 
as T — )■ 0; the modelled peak in Cm/T has the correct 
amplitude, however, it is broader and persists to higher 
temperatures. At low temperatures, the predicted inter- 
cept of Cm/T — ^ 7 decreases rapidly with increased de- 
fect density: the experimental intercept corresponds to a 
defect concentration of a; = 4.5% in our model. Since the 
linear specific heat arises from the long chain contribu- 
tions, alternatively, a modified chain length distribution 
([S]) that favours longer chains will match the experimen- 
tal intercept in Cm/T with a higher defect concentra- 
tion. This would correspond to a repulsion between de- 
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FIG. 3: Using the best fit parameters to the susceptibility 
[x = {5,6,8.5%}, J/™ = 34 K, Jafm = 63 K], the modelled 
specific heat of the broken chains and spin 1 defects. Note 
that the experimental intercept of Cm/T corresponds to a 
defect concentration of roughly x = 4.5% in our model. 



feet sites. The position of the peak and overaU scahng 
of Cm/T are determined primarily by the afm exchange 
in the Cu chains and would be largely unaffected by any 
correlations between defect positions. The measured heat 
capacity requires first a subtraction of the phonon contri- 
bution: its form was taken from the heat capacity of the 
non-magnetic isostructural compound Ba3ZnSb2 09 1]. 
Above T ~ 10 K, the phonon contribution overwhelms 
the magnetic response, and the result after subtraction 
may be unreliable; indeed, the experimental curve ex- 
trapolates rapidly to negative values in that regime (the 
dashed curve in Figure [s]). 

In summary, we propose that the magnetic degrees of 
freedom in Ba3CuSb209 decouple into antiferromagnetic 
spin chains due to a constant planar e^ orbital ordering. 
The defect Cu substitutions on Sb(2) sites couples to an 
in-plane Cu to form effective spin 1 orphans that simul- 
taneously break the spin chains at each defect site. Our 
modelled magnetic susceptibility reproduces the exper- 
imental curve remarkably well; however, the magnetic 
specific heat initial slope may require further consid- 
eration of correlations among the off-plane Cu defects. 
The predicted orbital ordering should to lead to an or- 
thorhombic distortion of the Cu-0 octahedral structure 
and should be discernible in further structural analyses 
of the compound. 
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